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ABSTRACT 


A  five  degree-of-freedom  linear  model  for  the  XR-3 
surface  effect  ship  is  developed  for  constant  speed  operation. 
The  weight  removal  transient  response  and  the  roll  behavior 
of  the  craft  in  turn  maneuvers  are  investigated  by  both 
linear  and  six  degree-of-freedom  nonlinear  models,  and  the 
linear  model  simulation  results  are  compared  with  the  non¬ 
linear  model  simulations.  Some  of  the  parameters  of  the 
craft  are  investigated  for  roll  motion. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  Surface  Effect  hybrid  type  high  performance  ship  has 
been  of  great  interest  recently,  and  in  1970  the  XR-3  test 
craft  was  delivered  to  the  Naval  Postgraduate  School.  As 
seen  in  17]  Leo  and  Boncal  converted  the  digital  simulation 
program  for  the  100-B  test  craft  to  represent  the  XR-3,  thus 
providing  a  nonlinear  model.  In  [1],  Gerba  and  Thaler  have 
developed  a  heave-only  linear  model  and  further  studies  were 
made  by  Barnes  in  [2]  by  developing  a  two  degree-of-freedom 
linear  pitch  and  heave  model  representing  the  XR-3  test  craft. 

B.  OBJECTIVE 

In  order  to  investigate  the  roll  behavior  of  the  XR-3 
test  craft  during  turn  maneuvers,  a  five  degree-of-freedom 
linear  model  is  developed  in  this  thesis.  Some  modifications 
are  made  for  pitch  and  heave  modes,  therefore  test  of  the 
vertical  plane  by  means  of  weight  removal  is  reinvestigated 
and  compared  with  the  nonlinear  simulation  results. 

In  [3] ,  the  roll  behavior  of  the  XR-3  and  the  sidewall 
sensitivity  are  studied  using  a  six  degree-of-freedom  non¬ 
linear  model.  Similar  studies  are  made  in  this  thesis  after 
developing  a  five  degree-of-freedom  linear  model  and  comparing 
with  the  nonlinear  model  results.  Finally,  some  modifications 
are  suggested  for  nonlinear  modeling. 

i 
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II.  NONLINEAR  XR-3  MODELING  FOR  FIVE  DEGREE-OF-FREEDOM 


A.  ASSUMPTIONS 

The  following  assumptions  and  simplifications  are  used 
during  development  of  a  constant  speed  XR-3  test  craft  model: 

1.  Aerodynamic  drag  on  the  forward  seal,  and  superstruc¬ 
ture  skin  friction  forces  are  neglected. 

2.  Momentum  (ram)  drag  is  disregarded  for  all  modes. 

In  other  words,  it  is  accepted  that  the  location  of  the  lift 
fans  with  respect  to  the  craft  C.G.  could  not  affect  signi¬ 
ficantly  the  roll  and  pitch  stability.  Also,  air  stagnation 
force  in  the  z-direction  is  assumed  to  be  negligible. 

3.  During  dynamic  response  of  the  craft,  only  air 
leakage  with  constant  leakage  area  at  the  stern  seal  is 
assumed.  No  other  air  escape  is  considered  such  as  possible 
leakage  from  sidewalls  or  bow  seal. 

4.  Plenum  pressure  and  air  mass  thermodynamics  are 
represented  by  ideal  diabatic  processes  with  the  air  constant 
specific  heat  ratio  y  =  1.4. 

5.  Effective  plenum  roof  area  is  developed  as  a  function 
of  craft  weight. 

6.  Planing  forces  are  concentrated  at  the  centroid  of 
the  keel  longitudinal  cross  section  through  the  fluid  medium. 

7.  Pitch  and  roll  damping  moments  due  to  added  mass  are 
included  in  the  craft  dynamics,  but  their  contribution  to 
heave  motion  is  neglected. 
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8.  The  thrust  forces  and  all  hydrodynamic  resistive 
forces  (drag)  in  the  2-direction  cancel  each  other,  since 

at  constant  speed  operation,  these  forces  are  in  equilibrium. 
However  they  might  contribute  to  the  moment  equations  such 
as  pitch  moment,  but  considering  their  effective  moment  arms 
with  respect  to  the  C.G.  are  almost  equal  to  the  thrust  moment 
arm,  so  that  their  contribution  in  the  moment  equations  is 
also  neglected. 

9.  For  seal  modeling,  vertical  resistive  forces  due  to 
hydrodynamic  response  are  taken  into  account  including 
corresponding  effective  wetted  seal  area. 

10.  Simplified  craft  geometry  is  used  for  plenum, 
buoyancy  and  seal  modeling. 

B.  •  COORDINATE  SYSTEM  AND  SIGN  CONVENTIONS 

1.  The  cartesian  coordinate  system  is  used  with  its 
origin  located  on  the  water  line  corresponding  to  the  craft 
center  of  gravity. 

2.  The  x-axis  lies  on  the  water  line  with  the  positive 
longitudinal  displacement  being  measured  forward. 

3.  The  y-axis  lies  in  the  vertical  cross  section  plane 
with  the  lateral  dislacement  being  measured  positive  to 
starboard. 

4.  The  z-axis  lies  in  the  vertical  cross  section  plane 
with  the  positive  vertical  displacement  being  measured 
downward. 
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5.  z  represents  heave  C.G.  displacement  measured  from 
the  origin  (at  water  line)  to  the  center  of  gravity. 

6.  z„  is  one  of  the  craft  dimensional  values  measured 

s 

from  the  keel  to  the  C.G. 

7.  Pitch  is  the  rotational  motion  about  the  y-axis, 
angular  pitch  displacement  is  represented  by  6 ,  positive 
pitch  angle  is  measured  upward  (boat  noses  up) . 

8.  Roll  is  the  rotational  motion  about  the  x-axis, 

<j>  angle  represents  angular  roll  displacement,  positive  roll 
angle  is  measured  to  starboard  (boat  heels  to  starboard) . 

9.  Similarly  the  yaw  is  the  rotational  motion  about  the 
z-axis,  positive  yaw  is  measured  in  the  clockwise  direction. 

10.  Positive  roll,  pitch  and  yaw  moments  are  to  be 
defined  as  the  moments  in  the  positive  roll,  pitch  and 
yaw  directions  correspondingly. 

11.  As  a  reference  zero  pitch  and  roll  angle  provide  the 
x-y  plane  parallel  to  the  water  surface. 

The  craft  coordinate  system  is  shown  in  figure  1. 

C.  SIMPLIFIED  CRAFT  MODELING  AND  GEOMETRY 

1.  General  Simplified  Geometry  and  Dimensions  of  the 
xr-3  Test  Craft 

In  figure  1  the  sideview  simplified  geometry  of  the 
XR-3  test  craft  is  shown.  Also  in  figure  3  the  top  view  is 
illustrated.  The  variables  and  dimensions  in  these  figures 
are  defined  as  follows; 
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geometry  ( top  view) 


I £  *  Draft  of  the  craft  at  center  of  gravity 

C  ■  Plenum  pressure  center  (the  centroid  of  the 
p  effective  plenum  area) 

*  The  longitudinal  distance  of  the  sections' 
centroids  from  the  C.G.  (forward  and  aft 
section  centroids  respectively) 

l~tl .  =*  Bow  and  stern  seals  front  face  distances 

at  water  level 

*  Centroid  distance  of  the  planing  force 
effective  area 

£.,*2  *  Forward  and  aft  sections  sidewall  submerged 

region  length. 

All  craft  dimensional  values  are  taken  from  the  XR-3 
six  degree-of-freedom  model  scale  drawings  as  indicated  in 
[2]  and  [3]. 

2.  Sidewalls 

As  reported  in  [2] ,  the  sidewall  modeling  is  developed 
from  the  SIDEWALL  subroutine  of  the  6  DOF  computer  model. 

Like  craft  dimensions,  all  sidewall  dimensions  and  their 
related  geometric  configurations  are  based  on  scale  model 
drawings.  For  simplicity  of  calculation,  the  sidewalls  are 
divided  into  two  sections  called  forward  and  aft.  During 
dynamic  response  of  the  sidewalls,  these  sections  are  con¬ 
sidered  individually  with  their  own  centroids,  so  related 
calculations  can  be  done  easily.  Especially,  the  sidewall 
modeling  is  developed  in  order  to  determine  buoyancy  forces 
and  moments.  By  means  of  this  modeling,  the  average  sidewall 
width,  average  draft  at  section  centroid  and  effective 
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plenum  area  are  also  determined  and  used  extensively  in  the 
rest  of  the  calculations. 

The  sidewall  vertical  cross  sections  at  the  section 
centroid  forward  and  aft  are  sketched  in  figures  4(a)  and  4(b) 
respectively.  From  these  figures,  the  average  forward  side- 
wall  width  at  the  buoyancy  center  is  found  to  be 


2  tan  6 . 


+  b  i 


(II-l) 


where 

=  Draft  at  forward  buoyancy  center 

tan  6.  *  Tangent  of  the  sidewall  deadrise  angle 

1  at  the  forward  buoyancy  center 

fa  =  Sidewall  flat  surface  width  at  the  kell 
f  (forward  section) 


Therefore  submerged  sidewall  volume  of  the  forward 
section  is  given  as  follows 


Ff  lf 


CII-2) 


can  be  formulated  as 

Af  *  Jld  -  *5  tan  8 

then  the  small  angle  approximation  tan  0*0  gives 
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Figure  5.  Plenum  modeling 
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I 


f 


(II-3) 


As  defined  before,  draft  at  C.G.  is 


la  ’  S  +  2 


( X 1—4 ) 


Similar  expressions  can  be  derived  for  the  aft 


section. 


2  tan  8-  +  ^a 


(II-5) 


where 


sidewall  flat  surface  width  at  the  keel 
(aft  section) 


*3  *2  ‘a 


(II-6) 


*d  +  *6  0 


Note  that  £_  are  functions  of  draft  and  pitch 

X  cl 

angle . 

3.  Plenum  Modeling 

In  [2]  plenum  modeling  is  presented  including  plenum 
pressure  variation  due  to  changing  of  draft.  In  the  develop¬ 
ment  of  a  pitch,  heave  and  roll  linear  model,  the  simplified 
plenum  geometry  and  dynamics  which  were  derived  by  Barnes 
are  considered  and  used  in  our  plenum  dynamics  representation 
The  variables  and  dimensions  for  plenum  modeling, 
as  seen  in  figure  5,  are  found  as 


2 


L»] 


( II— 8 ) 


£p  "  *8  “  Z9  (  h  } 


A .=  b  l 

O  W  p 


(II-9) 


=  ln  -  £ 

9  8  o 


where 

=  Length  of  plenum  at  the  water  line 

A^  =  Effective  plenum  area 

b  =  Lateral  distance  from  one  sidewall  to 
w  another  in  the  plenum  region. 

£g,  lQ  and  h  are  all  constants  related  to 

craft  dimensions.  Only  the  effective  plenum  area  A^  varies 
when  the  craft  draft  varies,  that  means  variation  of  the  plenum 
pressure  centroid  Cp  is  forward  or  backward. 

4.  Bow  and  Aft  Seals 

The  seals  are  an  involved  additional  dynamic  system 
with  their  own  degree  of  freedom  relative  to  the  craft,  and 
they  carry  little  load  compared  with  the  other  components 
of  the  craft.  However  in  the  simplified  model,  the  seal 
motions  relative  to  the  craft  are  neglected. 

The  assumed  seal  geometry  is  shown  in  figure  2, 
developed  in  [2].  Using  an  appropriate  geometric  configura¬ 
tion,  the  seal  wetted  area  is  derived  as  a  function  of  craft 
draft  and  pitch  angle. 
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From  [2],  the  width  of  the  wetted  surface  for  the 


bow  seal  is; 


£.  -  £,  tan  8 
xf  "  - sin~~5'P - 


(11-10) 


where 


sin  31°  =  sin  of  the  angle  between  water  surface 

and  the  bow  seal  front  face. 


The  numerator  of  this  ratio  represents  draft  at  the 
bow  seal  front  face. 

A  similar  expression  for  the  stern  seal  is; 


x 


a 


+£^  tan  8 
sin  32° 


Multiplying  xf  and  xa  by  fa w,  bow  and  stern  seals 
wetted  areas  are  obtained  as  follows; 


=  fa. 


w 


=  b 


w 


5.  Rudder  and  Propulsion  Hydrodynamics 

The  rudder  hydrodynamics  is  calculated  from  [4] , 
inserting  rudder  geometry  dimensional  values  into  the  fin 
equations.  Actually  we  are  interested  in  the  lateral  force 
(sway)  on  the  rudders  for  a  given  rudder  order  deflection. 
This  dynamic  equation  is  given  in  [4]  as  follows; 
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(Il-ll(a) 


-  »  “2  Ar(1  +K^) 


d  Vh 

crl6r-(ltE;>  V1 


where 

Ar  *  Rudder  platform  area  (for  one  rudder  only) 

dg  =  draft  at  the  stern  =  z  +  zg  -  xr  0 

hr  =  Vertical  depth  to  bottom  tip  of  the  rudder 

Cr  =  Rudder  lift  coefficient 
2  it  R 

=  ^ j  (for  fully  wetted  flow  regime) 

d 

R  *  The  rudder  aspect  ratio 

6  =  Rudder  deflection  (the  sign  convention  of  the 

r  rudder  is  to  be  positive  for  port  turn) 

vh  *  V  +  xr  r  -  2S  p 

x  =  longitudinal  coordinate  of  rudder  geometric 
r  center  (always  negative  value) 

V  =  lateral  speed  (ft/sec) 

r  *  yaw  rate  (rad/sec) 


In  order  to  simplify  the  rudder  dynamic  equation, 
the  following  coefficients  are  defined  as 


2  ds 

ct  -  P  u  Ar(l  +  ji)  Cr 
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Then  equation  (Il-ll(a))  becomes 


Y 


r 


Ct(6r 


(Il-ll(b)l 


Let  us  define  zr  to  be  the  vertical  coordinate  of 
the  rudder  from  the  C.G.  to  its  geometric  center,  the  rudder 
roll  and  yaw  moments  can  be  calculated  from 


Kr  =  -  Yr  zr  (11-11 (c) ) 

Nr  =  Yr  xr  (.11-11  (d) ) 


During  turn  maneuvers,  the  propulsion  thrust  vector 
is  also  deflected  according  to  the  given  rudder  order.  There¬ 
fore  the  lateral  component  of  the  thrust  contributes  yaw  and 
roll  moments  as  well  as  a  sway  mode,  but  this  diturbance  must 
be  considered  as  a  step  input  to  the  linear  modes  with  constant 
magnitude,  therefore  the  total  propulsion  sway  force  is 


Ypr  =  2  Xfc  sin  5r  (11-12) 

where  is  the  propulsion  thrust  corresponding  to  craft  speed. 

Similarly  the  roll  and  yaw  moments  generated  by  the 
propulsion  system  are 


K 


Pr 


(11-13 (a) ) 
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(11-13 (b) ) 


where  xp  and  zp  are  the  longitudinal  and  vertical  locations 
of  the  propeller  respectively. 

D.  DYNAMIC  FORCES  FOR  FIVE  DEGREE-OF-FREEDOM  MODES 
1.  Seal  Forces 

Considering  the  seal  configuration  and  dynamics,  the 
pressure  difference  across  the  wetted  portion  of  the  seal  is 


where 

_  2 
Pb  =  Plenum  gauge  pressure  (lb/ft  ) 

Pa  *  Atmospheric  pressure 


Since  the  seal  is  in  equilibrium,  the  vertical  hydro- 


where 
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where 


A  =  flat  surface  area  at  the  keel 
u  =  speed  of  the  craft 

Using  the  small  angle  approximation  and  including  the 
loss  coefficient  Cpl  suggested  by  Barnes  in  [2] ,  the  planing 
force  is  redetermined  as 

Zpl  =  -  Cpl  it  p  A  u2  0  (11-18) 

A  similar  expression  is  derived  for  sidewall  dynamics 
in  [4],  but  in  our  modelling  equation  (11-18)  is  used. 

4 .  Pitch  Damping  Forces 

The  sidewalls  of  the  surface  effect  ship  influence 
moderately  heave,  pitch  and  roll  modes,  so  in  (4]  the  hydro- 
dynamic  and  hydrostatic  forces  and  moments  acting  on  the  side- 
walls  of  uniform  triangular  cross  section  are  derived  with 
particular  application  to  vertical  plane  motion  and  lateral 
motion. 

Considering  constant  speed  operation  of  the  craft, 
neglecting  some  terms  in  the  equations  provided  from  [4], 
the  vertical  force  acting  on  only  one  sidewall  is  found  as; 

z,w1  ■  A33  xs  u  *  (ii-19> 
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where 


=  vertical  added  mass  calculated  at  the  stern 
=  (ir/8)  p  Bs2 

B  =  value  of  sidewall  waterline  at  the  sidewall 
stern 

xg  =  stern  distance  from  C.G.  (negative  value) 
q  =  pitch  rate  (rad/sec). 

For  the  pitch  damping  moment,  the  same  equation  is 
used  by  Barnes  in  [2]. 

5.  Roll  Damping  Forces 

The  XR-3  test  craft  is  dynamically  similar  to  100-B 
surface  effect  ship.  Therefore  we  can  use  dynamic  equations 
calculated  in  [4]  for  the  XR-3  test  craft.  The  heave-roll 
damping  force  for  one  sidewall  can  be  written  as; 

zsw2  *  -  *33  u  bo  P  (II-20) 

In  the  above  equation,  faQ  is  half  of  the  craft  width 
and  it  must  be  used  (+)  for  starboard  sidewall,  (-)  for  port 
sidewall. 

Similarly  the  lateral  damping  force  due  to  sidewall 
hydrodynamics  is 


YSWl  =  "  A22  u(v  +  xsr  ‘  zsp)  (II_21) 
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here  A22  is  the  lateral  added  mass  calculated  at  the  stern, 
p  is  the  roll  rate.  A22  is  calculated  from  [4] ,  that  is 

2 

A22  =  tt/5  p  dg 

6.  Plenum  Pressure  Lift  Force 

From  [2]  and  [4],  the  plenum  pressure  lift  force  is 
determined  assuming  uniform  pressure  distribution.  This  is 
formulated  as 


5p  * 


from  equations (II-8)  and  (II-9)  follow 


Zp  *  -  fcw  V*S  - 


(11-22) 


E.  DYNAMIC  PLENUM  AIR  MASS,  PRESSURE  EQUATIONS 

Assuming  constant  aft  seal  leakage  area,  the  leakage  flow 
rate  (Qout)  is  given  by  Gerba  and  Thaler  for  the  XR-3  heave 
only  model  in  [1] ,  then  we  have 


2  P 


‘out 


-  a,  ( 


n  "4'p. 


b)  1/2 


(11-24) 


where 

C  ®  Nozzle  flow  coefficinet 
n 

Ag  =  Leakage  area 
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Atmospheric  air  density. 


f 


The  linearized  air  flow  (Q^n)  into  the  plenum  is  pre¬ 
sented  in  [5]  according  to  the  fan  properties  and  calculated 
in  (1]  as  follows. 


‘in 


n(Q0  " 


V 


(11-24) 


where  n  is  the  number  of  fans  supplying  the  air  to  the  plenum, 
Qo  is  the  steady  state  air  flow  rate  corresponding  to  steady 
gauge  pressure,  is  (1/slope)  of  the  vs.  Q^n  curve  at 
the  equilibrium  point. 

From  equations (11-23)  and  (11-24),  the  plenum  air  flow 
rate  can  be  written  as 


“b 


5b  _ 

dt 


pa ^®in  ~  ®out* 


(11-25) 


The  equation  (11-25)  reveals  thatif  the  fan  flow  rate 
(Qin)  equals  to  the  leakage  flow  rate  (Qout) »  then  the  plenum 
air  mass  (M^)  is  constant. 

The  absolute  plenum  pressure  is  calculated  from 


-  Vs 


«b 


(11-26) 


In  the  above  equation,  is  defined  as  plenum  volume 
and  its  simplified  form  is  calculated  from  [1J , 
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<vn  -  *b  V 


(11-27) 


where  Vn  is  to  be  the  empty  plenum  volume. 

F.  NET  HEAVE  FORCE 

The  resultant  heave  force  is  the  summation  of  the  all 
external  dynamic  forces/  it  is 


Z  -  +  Z  +  Z.,  +  Z,  +Z.+Z  +  Z 

sf  sa  bf  ba  pi  sw  p 


where 


Z  *  Z  +  Z 

SW  SWj^  sw2 

Zgw  also  includes  a  pair  of  heave  buoyancy  forces  due  to 
roll,  but  they  are  in  opposite  directions.  Therefore  these 
stabilizing  forces  are  only  taken  into  account  in  the  roll 
motion. 

G.  PITCH  AND  ROLL  MOMENTS 
1.  Seal  Moments 
a.  Pitch  Mode 

The  seal  force  is  modeled  by  Barnes  in  [2]  acting 
on  the  center  of  the  seal  wetted  length.  The  effective  pitch 
moment  arms  are; 

£gf  *  £3  -  (for  forward  seal)  (11-28) 

£  -  l.  +  x./2  (for  aft  seal)  (11-29) 

sa  4  a 
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From  this  follows  the  forward  and  aft  seal  moments 


M 


sf 


Z  -  i  , 
sf  sf 


M 


sa 


sa  sa 


(11-30) 

(11-31) 


where  Z  _  and  Z  have  already  been  defined, 
b.  Roll  Mode 

The  flexible  seals  can  contribute  little  to  roll 
stability.  Actually  the  seal  configurations  and  hydrodynamics 
are  relatively  more  complex.  In  (4],  the  seal  roll  moments 
are  calculated  by  combining  the  roll  contribution  of  all  seal 
elements. 

Considering  the  seal  response,  the  submerged 
region  of  the  seal  can  be  calculated  and  its  roll  moment  con¬ 
tribution  can  be  formulated  in  the  following  equations; 

Ksf  -  -  6W3/12  p  g  xf  tan  ♦  di-30) 

Ksa  =  “  6W3/12  p  9  xa  tan  *  (11-31) 

For  this  modeling,  the  roll  moment  arm  is  considered 
as  (2/3  faQ) ,  that  is  the  centroid  of  the  triangle  cross  section 
of  the  submerged  portion. 

2.  Buoyant  Moments 
a.  Pitch  Mode 

The  pitch  buoyancy  moments  are  calculated  from 
equation  (11-16)  and  (11-17)  assuming  constant  moment  arms 

3  2 


&5  and  £g 


The  pitch  buoyancy  moments  are 


Mtf  -  -  Zbf  t5  (11-32) 

«ba  ’  Zba  “5  'IJ-33> 


b.  Roll  Mode 

As  discussed  in  the  heave  mode,  the  effect  of  the 
sidewall  buoyancy  due  to  the  roll  is  to  be  a  stabilizing  moment. 
From  figures  4(a)  and  4(b)  the  buoyancy  volume  change  can 
be  approximated  as 


AV^  = 


bf+bw 
(■■?--*  + 


2  tan  6 


-)  ( 


tan  B;  +  bf)4l* 


(11-34) 


Finally,  roll  moment  is 


Kbf  =  -  2  p  g  AVf  (11-35) 

Similarly,  the  aft  section  roll  moment  is 


*ba 


-  2  p  g  AVa 


3.  Planing  Pitch  Moment 

The  planing  effect  due  to  the  sidewall  flat  area  at 
the  keel  contributes  only  pitch  motion  as  follows; 
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4.  Pitch  Damping  Moment 

This  is  calculated  from  equation  (11-19)  as 


M 


sw 


-  2  Z 


sw. 


-  2 


A_ ,  u  x 
33  s 


In  the  above  equation,  the  factor  "2"  appears  because 
there  are  two  sidewall  sections. 

5.  Roll  Damping  Moment 

Combining  equation  (11-20)  and  (11-21),  the  roll 
damping  component  is 


'sw  *  2  2sw2  bo  -  *sw1  2s  *  “wave  (II-37> 


where  Kwave  represents  the  additional  roll  damping  tern 
because  of  vertical  wave  generation  in  the  roll. 

From  [3]  and  [4]  ,  Kwave  is 


wave 


=  -  32/tt  b 


Bc  P 


(11-38) 


Here,  Bc  is  the  sidewall  vertical  added  mass  given  as 

2  2 

IS  =*  -rr/8  p  l  i .  Cotan  8  . 

c  a 


6.  Plenum  Pitch  Moment 

From  equation  (11-22)  and  defining  Xcp  to  be  the 
plenum  pressure  centroid,  the  plenum  pitch  moment  is  given 
as 

i 

i 
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«p  -  -  2p  Xcp  "  VPb  -  Pa>  xc 


where  X  is  calculated  as  a  function  of  craft  draft  in  [2] 
cp 

as  follows; 


X  =  X  +  (i  i  ,)/2h 
cp  c  9  d  ' 


X  =  IJ 2  +  K 

c  o  c .  9  • 


7.  Plenum  Roll  Moment 

The  destabilizing  static  plenum  roll  moment  is 
generated  when  the  craft  rolls  to  one  side.  Since  the  effec¬ 
tive  plenum  pressure  centroid  shifts  laterally  to  the  other 
side  tending  to  increase  rolling  of  the  craft.  The  plenum 
roll  moment  is 


'p  ■  Vpb  *  Pa)Us  '  V* 


8.  Rudder  and  Propulsion  Roll  Moments 

These  roll  moments  have  already  been  calculated 
from  equation  (II-ll(c))  and  (11-13  (a)). 


H.  SWAY  FORCES 


1 .  Cross  Flow  Drag  (Sidewall) 


During  the  turn  maneuver,  the  cross  flow  drag  is 
generated  acting  on  the  outboard  sidewall.  The  direction 
of  this  force  is  negative  while  turning  to  the  left. 
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Neglecting  the  lateral  drag  component  inside  of  the  sidewall, 
the  net  cross  flow  drag  can  be  calculated  from 

Dsw  =  '  I  p  si9n<V  Cd  Vr2  M*d  4  fao4>)  (XI"39) 

where 

vr  =  V  +  xs  r  -  Cig  -  V 2)P 

=  cross  flow  drag  coefficient 

Considering  the  outboard  surface  of  the  submerged 
sidewall  sections  (dead  rise  angle  effect) > the  lift  force 
in  the  z-direction  can  be  produced  due  to  the  lateral  drag 
deadrise  projection.  In  [4]  ,  the  importance  of  the  sidewall 
deadrise  angle  is  mentioned,  however  at  the  hydrodynamic 
section  in  this  reference,  the  lift  component  of  the  drag 
force  has  not  been  derived. 

In  this  modeling,  the  lift  force  (Zd)  and  lateral 
drag  component  is  calculated  assuming  both  have  the  same  drag 
coefficients  (Cd  =  1.28).  They  are 

Y.  *  D  sin  0  (11-40) 

d  sw 

Z.  »  Doi_r  cos  6  (11-41) 

a  sw 
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From  equation  (11-40)  and  (11-21),  the  total  hydro 
dynamic  sway  force  of  the  sidewall  is 


sw 


=  (D  +  Y  )  sin  6  +  Y 


sw 


sw. 


sw. 


(11-42) 


2.  Centrifugal  Force 

The  centrifugal  force  acting  on  the  C.G.  of  the  craft 
is  always  in  the  outboard  direction  during  the  turn  maneuver. 
Yc  is  the  centrifugal  force,  given  as 


Y 


c 


R,  the  turning  radius  is  equal  to  (u/r) ,  combining  this 
with  the  above  equation  gives 


Y_  =  -  m  u  r 

c 


(11-43) 


a  negative  sign  is  used  due  to  the  coordinate  system  sign 

convention.  Since  the  right  turn  produces  (+)  yaw  rate,  then 

Y  is  calculated  as  (-) ,  so  outboard  sway  force  is  generated 
c 

for  this  configuration. 

3.  Rudder  and  Propulsion  Sway  Forces 

They  are  determined  from  equation  (Il-ll(b))  and 

(11-12) . 

The  sway  forces  due  to  the  seals ,  air  plenum  and 
aerodynamic  drag  are  neglected  in  this  modeling. 
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I .  YAW  MOMENTS 


1.  Sidewall  Yaw  Moment 

The  sidewall  yaw  dynamic  equations  are  given  in  [4]. 
After  some  simplification,  the  net  widewall  yaw  moment  is 
calculated  from 


N 


sw 


N1  +  N2  + 


(X  ,  -  X  )  b 

pw  sw  o 


(11-44  (a) ) 


where 

=  2u(v  + xgr  -  zgp)A22£ (1  -  X) 

-  tt/5*  p  (Jld  +  bo4> )  2£u  (v  -  zsp)  (11-44  (b)  ) 

N2  =  -  p  Cd  £d(2  S1  £3  r  v  +  s2  l2  Vr2)  (11-44 (c) ) 


Z  +  X 

X  =  — - — =  The  ratio  of  C.G.  location  of  (11-44  (d) ) 

sidewall  bow  to  the  craft 
length 


_  1  -3M1  -  X) 

l  ~  3 


(11-44  (e) ) 


s2 


(11-44 (f ) ) 


-  1/2  p  C.  A  u2  (11-44 (g)) 

P 

port  wetted  area  surge  resistive  force 


X 


sw 


-  1/2  p  Cf  Aw^  u2 


(11-44 (h) ) 


=  starbord  wetted  area  surge  resistive 
force 
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=  Average  sidewall  flat  surface  width  at  the  keel. 

In  equation  (11-45)  ±  term  related  to  the  starboard 
or  port  wetted  surface  calculations. 

2.  Rudder  and  Propulsion  Yaw  Moments 

The  yaw  moments  due  to  the  rudder  and  propulsion  sys¬ 
tem  are  determined  from  equation  (Il-ll(d))  and  (11-13 (b) ) 
respectively. 

The  seals,  air  plenum  and  aerodynamic  yaw  moments  are 
neglected.  Actually  the  aerodynamic  yaw  moment  should  be  taken 
into  account.  Since  during  the  turn  maneuver,  the  yaw  drift 
angle  is  produced,  thus  a  positive  aerodynamic  yaw  moment 
contribution  to  the  yaw  motion  has  been  generated.  From  6 
DOF  nonlinear  model  simulations,  this  contribution  is  determined 
at  5.2%  of  the  total  yaw  moment. 

J.  TOTAL  PITCH  MOMENT 

The  summation  of  the  all  dynamic  pitch  moments  gives  the 
net  pitch  moment  which  can  be  written  as 


K.  TOTAL  ROLL  MOMENT 

Before  calculating  net  roll  moment,  the  roll  moment 
equation  due  to  the  cross  flow  drag  term  must  be  derived. 
Including  also  the  added  mass  effect  to  the  equation  (11-40) 
and  (11-41) ,  the  result  is 

Kd  =  (Dsw  +  Ysw1) (fao  cos  B  ’  zo  Sin  S)  (11-47) 

Combining  all  roll  dynamic  equations,  the  net  roll 
moment  is  calculated  as 


K  =  K 


sf 


+  K  + 
sa 


Kbf  + 


K,  _  +  K  +  K  +  K  + 
ba  sw  pr  r 


K  +  K, 
P  a 


L.  TOTAL  SWAY  FORCE 

From  equation  (Il-ll(a)),  (11-12),  (11-42)  and  (11-43), 
the  total  sway  force  is  obtained  as 


Y  = 


Y  +  Y  +  Y  +  Y 
r  pr  sw  c 


(11-48) 


M.  TOTAL  YAW  MOMENT 

The  summation  of  the  sidewall,  rudder  and  propulsion 
yaw  moments  gives  the  net  turning  moment  about  the  z-axis 


N 


N  + 
sw 


Nr  + 


N_ 


Pr 


(11-49) 
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N.  FIVE  DEGREE-OF-FREEDOM  EQUATIONS  OF  MOTION 


Considering  equations  of  motion  for  XR-3  test  craft, 
the  yaw  is  defined  as  the  rotational  motion  of  the  craft 
coordinate  system  relative  to  the  fixed  frame  and  the  positive 
\p  angle  is  measured  in  the  clockwise  direction.  This  sign 
convention  is  in  agreement  with  the  rudder  order  sign  defini¬ 
tion.  Since  positive  rudder  angle  corresponds  to  turning 
to  the  left,  producing  negative  yaw  moment  and  finally  at 
the  equilibrium  condition,  the  negative  yaw  moment  is 
generated. 

Assuming  small  values  of  pitch  (0)  and  roll  angles  and 
relating  heave  motion  to  the  free  surface,  the  relationship 
between  the  fixed  and  craft  coordinate  system  can  be  written. 
Eliminating  the  terms  related  to  the  fixed  coordinate  system, 
the  simplified  5  DOF  equations  of  motion  are  obtained  as 
follows; 


m  v  =  Y 


m  w  =  Z 


I 


XX 


p 


I 

xz 


r 


K 


I 


zz 


r  - 


xz 


N 


(11-50) 

(11-51) 

(11-52) 

(11-53) 

(11-54) 
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The  right  hand  sides  of  the  equations  of  motion  give  the 
forces  and  moments  acting  on  the  craft  in  terms  of  the  craft 


in  terms  of  the  craft  coordinate  system  variables.  In 
equation  (11-57)  w  represents  heave  acceleration  and  all 
other  variables  in  these  equations  have  already  been  defined. 

In  order  to  develop  a  5  DOF  linear  model,  only  the  surge 
equation  of  motion  (x-direction)  is  neglected.  Since  from 
6  DOF  nonlinear  simulation  results,  very  small  velocity 
reduction  is  observed  during  the  turn  maneuver  of  the  XR-3 
test  craft  (over  the  speed  range  20  through  30  knots). 


III.  LINEARIZATION  OF  THE  EQUATIONS 

A.  TAYLOR  SERIES  EXPANSION  AND  TOTAL  DIFFERENTIAL 

Generally,  the  linearization  process  can  be  performed  by 
means  of  a  Taylor  series  expansion  using  the  linear  part  of 
the  series  (omitting  high  order  terms).  The  same  process  can 
be  easily  done  applying  a  total  differential  to  the  net  force 
and  moment  equations,  since  the  effect  of  small  disturbances 
to  the  system  is  to  be  determined. 

The  total  differential  of  a  given  dependent  variable  A 
with  respect  to  the  independent  variables  x,  y  and  z  is  given 
as 


-  f£  d*  +  If  di-  +  H  dz 

B.  LINEARIZATION  OF  THE  EQUATIONS  OF  MOTION  AND  AIR  FLOW 
DYNAMICS 

The  total  differentials  of  the  equations  of  motion  and 
air  flow  dynamics  have  been  derived  as 


dw  =  dZ/m 
dq  =  dM/Iyy 


dm.  -(pa  n  k.  +  C  A.  ( - 2__)1/2JdP. 

*  a  P  n  £  2pb(Q)  t 


dv 


dY/m 


(III-2) 

(III-3) 

(III-4) 

(III-5) 
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where  Pb(0)  represents  the  steady  state  value  of  the  plenum 

gauge  pressure  according  to  the  final  equilibrium  condition. 

The  equations  (11-52)  and  (11-54)  both  involved  second 

order  yaw  and  roll  acceleration  terms.  Therefore  applying 

the  total  differential  and  solving  these  two  equations 
■  • 

for  dp  and  dr  we  have 


dr  =  dK/I  +  dN/Ir  (III-6) 

dp  =  dK/Ik  +  dN/I  (III-7) 


I,  Ir  and  1^  are  defined  as 


A  1  1  "  1 

I  =  XX  22  X2 


XZ 


I  I  -  I 

A  XX  ZZ  XZ 

*xx 


A  1  I  -  I  ‘ 

A  XX  ZZ  XZ 

-  - J - 

ZZ 


(III-8  (a) ) 


(III-8 (b) ) 


(III-8  (c) ) 


Equation  (III-6)  and  equation  (III-7)  contain  both  net 
incremental  roll  (dK)  and  yaw  (dN)  moments.  This  is  simply 
a  cross  product  inertial  term  (I  )  effect  in  the  equations 
of  motion. 

In  equation  (11-59),  dP^  is  not  a  state  variable,  it  can 
be  considered  as  an  output  variable,  therefore  from  equation 
(11-26)  and  (11-27),  dPb  must  be 
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3PK  3Ph 

dPb  =  TT  dz  +  ^  ^ 

C.  FIVE  DEGREE-OF-FREEDOM  LINEAR  MODEL  SENSITIVITY  COEFFICIENTS 

The  sensitivity  coefficients  related  to  the  right  hand  side 
of  the  equations  of  motion  are  defined  as  the  partial  deriva¬ 
tives  of  the  net  forces  and  moments  divided  by  craft  mass  or 
inertial  terms. 

In  [2],  Barnes  has  developed  a  2  DOF  (heave  and  pitch) 
linear  XR-3  model.  Therefore  most  of  the  neave  and  pitch 
sensitivity  coefficients  are  directly  used  in  order  to  develop 
a  5  DOF  linear  model.  However  some  modification  is  made  for 
plenum  sensitivity  coefficients  and  pitch  mode.  Also  some 
of  the  steady  state  operating  point  values  have  been  reinves¬ 
tigated  using  the  6  DOF  nonlinear  model. 

In  all  the  linearized  equations,  the  final  steady  state 
values  of  the  variables  have  been  used. 

1.  Pitch  and  Heave  Sensitivity  Coefficients 

In  [2]  equation  (III-40) ,  the  sensitivity  coefficient 
due  to  z  has  been  redefined  as 


DHPZ 


b  .  cr 


-  V°>  +  p: 


This  is  claculated  from  (11-22)  where  the  second  term 
does  not  appear  in  [2]. 

The  modified  form  of  the  plenum  pitch  sensitivity 
coefficient  due  to  state  variable  z  is  to  be 
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OPPBZ 


^9  —  ^  ^9 

*  -  V0)  2H  -  xcp(0)  V0)  -nr1 

The  derivation  of  these  modified  equations  is  based 
upon  the  total  differential  of  the  plenum  lift  force,  it 
is 

3Z  _  3Z 

dZ  =  — E-  dP,  +  dA. 

p  3Pb  b  “b  ^ 

From  equations  ( IX— 9 ) ,  (11-26)  and  (11-27),  the  total 
differentials  dVfa  and  dPfa  have  been  calculated  as 

faw  £9 

dVb  =  ~(Ab(0)  -  Zd(0)  — j~)  dz  (III-9) 

®b  *  i'pb(0»I§or  +  v^Tor(Ab<«-*d(0>  (III-l0> 

dPfa  and  dVfa  equations  are  also  different  from  the  equations 
derived  in  [2]. 

The  bow  and  aft  seal  pitch  moments  sensitivity 
coefficients  (DPSFTH  and  DPSATH)  due  to  0  and  plenum  air 
mass  equation  (AMB)  are  rederived  and  their  modified  forms 
are  used  for  the  pitch  mode  in  the  5  DOF  linear  model  computer 
program. 

Considering  the  5  DOF  linear  model  state  space  matrix 
representation,  the  a^'s  represent  "A"  matrix  elements. 

Therefore  the  a^j  values  relating  to  the  sensitivity  coeffi¬ 
cients  for  heave  motion  are  determined  as 
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DZZ 


a21  = 
a23  =  DZTH 

a24  =  DZDP 

a25  =  DZMB 

where  DZZ,  DZTH,  etc.,  have  been  derived  in  [2].  However 
their  modified  forms  can  be  found  in  the  5  DOF  linear  computer 
program. 

For  pitch  motion,  a^j  coefficients  are 

a..  =  DTHZ 

41 

a..  =  DTHTH 

43 

a..  »  DTHDTH 

44 

a.  c  =  DTHMB 

45 

2.  Plenum  Air  Flow  Sensitivity  Coefficients 

These  coefficients  are  calculated  directly  from  [2]. 
The  a^j's  related  to  plenum  air  sensitivity  coefficients  are 

a51  =  DMBZ 
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DMBMB 


a55  " 

3.  Sway  Sensitivity  Coefficients 

Before  deriving  the  sensitivity  coefficients  for 
sway  motion,  some  partial  derivatives  related  to  the  rudder 
and  sidewall  hydrodynamics  must  be  calculated.  From  equation 
(11-11 (b) )  Yr  is  already  known,  this  follows  the  partial 
derivatives  of  Yr  with  respect  to  the  state  variables 

=  Ct<Ce/u)  Zg  (III-ll) 

=  -  Ct(Ce/u)  (III-12 ) 

=  "  C.  (C  /u)  X  (III-13 ) 

t  e  r 

From  equation  (11-39),  the  sidewall  cross  flow  drag 
terms  D  partials  are  derived  in  the  following  equations. 

3W 

Let  us  define  first  some  coefficients 

di  4  r»cd‘ 

dr  "  +  si9n(Vr)  faQ  ( 0 ) 

The  partials  are 

3D  - 

Dz  =  =  ‘  di  sign(Vr)  v/(0)  (III-14) 


R_  = 


3Y 

_ i 

3p 


R 


3Y 

r 


R 


3Y 

r 

Tr~ 
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FT— 

D*  -  -  dl  Vr2<0>  bo 

(III-15) 

Dp  =  2  d1  dr  sign(Vr)  Vy (0)  (Zg  -  £d(0>/2) 

(III-16) 

Dv  =  "  2  dl  dr  V0) 

(III-17) 

Dr  =  -  2  di  dr  sign(Vr)  Vr(0)  Xg 

(III-18) 

d2  and  d^  are  defined  as 

d2  =  2  A33  bo  u 

(III-19 ) 

A  . 

d3  A22  U 

(III-20 ) 

J 

where  u  is  the  craft  speed  assumed  to  be  constant. 

For  sway  mode,  the  state  matrix  elements  are 

1  AY  1 

*96  =  i<!?>  =  £(D*  Sin  *> 

(III-21) 

where  is  given  from  equation  (III-15) . 

i .  3Y. 
a97  ml9p' 

=  |[R  +  (R  +  d,Z  ) sin  6  +  d,Zl 

m  p  p  -3  s  j  s 

(III-22) 
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_ A 

r 


i  ,ay. 

a99  '  m^v' 

=  ^[Rv  +  (Dv-d3)sin  6  -  d3J  (III-23) 

a  = 

911  m  3r 

=  ^[Rr  +  (Dr  -  d3xg)  sin  g  -  d3Xg -mu] 

( III-24 ) 

4 .  Roll  and  Yaw  Sensitivity  Coefficients 

Recalling  equations  (III-6)  and  (III-7) ,  the  net 
roll  and  yaw  moments  appear  in  both  equations  of  motion.  In 
this  case,  to  determine  sensitivity  coefficients  is  not  easy, 
since  the  yaw  and  roll  partials  are  involved  in  the  dynamic 
equation  with  different  inertial  terms  and  coefficients. 

In  order  to  simplify  the  problem,  the  roll  and  yaw 
partial  derivatives  are  calculated  individually,  then  combining 
these  coefficients  and  dividing  by  corresponding  inertial 
terms  give  the  system  sensitivity  coefficients  for  the  roll 
and  yaw  modes. 

a.  Roll  Mode 

The  net  roll  partial  derivatives  are  defined  by 
variables. 

The  partial  derivative  with  respect  to  z  is 

k71  -  H  ’  &(Ksf  +  Ks*  +  Kbf  *  "W  (III-25> 

In  equation  (III-26)  all  roll  components  have 
already  been  derived  in  the  preceding  sections . 
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From  equations  (11-30)  and  (11-31) ,  the  seal 


partials  are 


-  (573) 3  P  g  (1/sin  31°)  <»(0) 


-  (273) 3  P  g (1/sin  32°)  0(0) 


For  buoyancy  partials,  let  us  define  the  following 


coefficients ; 


b,  =  2  p  g  £. 


(III-26) 


b0  =  2  p  g  i. 


(III-27) 


.  4  ,6£  *  V  .  *f 

b-j  -  ( - >3 - )  +  o — T~r 


2  tan  6- 


(III-28) 


b4  k  tanTJ  +  faf 


(III-29) 


bc  4  ,-JLa—  +  fa 

6  tan  02  a 


(III-30) 


From  equations  (11-35)  and  (11-36) ,  the  buoyancy 


partials  are 


3Kbf 


bl  b3 

tSH-§7(b5+b3>  *<0) 


b2  b4 

ti5-^(b6+b4>*<01 
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Equation  (11-26)  follows; 

3K  _  3P. 

s?  -  -1^(0)  Pb(0)  +  (zs-Jd(0)  5/]  *(0) 

Similarly  for  sidewall  cross  flow  drag  partial  is 

3K, 

JT  =  Vfao  cos  3  '  Zo  sin  3) 

The  roll  partial  derivatives  with  respect  to 

0  are 


where 


k73 


3K 

30 


3T(Ksf  +  Ksa  +  Kbf  +  Kba> 


(III-31) 


3K 


sf 


30 


(2T3 ) 3  p  9<£3/sin  31°)  *(0> 


3Ksa 

TT 


-  (27j)3  p  g(Ji4/sin  32°) 


9Kbf 

T5~ 


tS5T7(b5  +  b3  V  *(0) 


5Kb« 

30 


b2  b4 


(b,  +  b.  i,)  4>(0) 


tan  32'“6  ~4  ~3‘ 
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is 


The  roll  partial  derivative  with  respect  to 


'75 


9K 


9P. 


=  ^  (0)  (Zg  -  Jld(0))  ^  <)>(0) 


(III-32 ) 


The  roll  partials  with  respect  to  <p  are 


where 


9K 

76  90 

=  -^r(K  -  +  K  +  K,  ,  +  K,  +  K  +  K,) 

90'  sf  sa  bf  ba  p  d 

(III-33) 


3Kbf 


3K 


ba 


9$ 


-  (bw/2.3)3  p  g  Xf 

-  (by 2. 3)  3  p  g  Xa 

-  bl  b32  b5 

-  b2  b42  b6 

Pb(0)  (Zs  -  £.d«0)  )  l 
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3K. 

JT  =  V6o  cos  6  "  Zo  sin  8) 


The  roll  partials  with  respect  to  p  are 


3K 

77  3p 

=  4-(K  +  K .  +  K  )  CHI-34) 

3p  sw  a  r 


where 


3K 

sw 

~W 


-  d_  -  Z  Z  Ao0  u  -  32/n  b  Z  B 
4  S  2  ^  o  c 


y 

3K  sw, 

3p~  "  (Dp  +  Ip— >  (bo  cos  8  '  zo  sin  81 


where 


The  roll  partials  with  respect  to  V  are 


L79 


3K 

3V 


W  (Ksw  +  Kd  +  V 


(III-35 ) 


3Kd 

W~  ~ 


(Dv  -  d3)  ( bQ  cos  3 


Z  sin  3)  +  d_  X  Z 

O  3  S  S 


3K 

sw 

3V 


d3  Zs 
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V* 


W~  =  ct(Ce/u)  Zr 


The  roll  partials  with  respect  to  yaw  rate  are 


k711  " 

3K 

ar 

= 

£(Ksw  +  Kd  +  Kr) 

(III-36) 

where 

9Ksw  _ 
3r 

d3  Xs  Zs 

3K 

sw 

ar 

(Dr  ‘  d3Xs> (6o  °°s 

5  -  ZQ  sin  6)  < 

3Kr 

3r 

Ct(Ce/u)  Xr  3r 

b. 

Yaw  Mode 

The  following  coefficients  are 

defined  for  yaw 

rotational  i 

motion; 

nl 

=  2  A 

A22 

Ml  -  A) 

(III-37) 

n2 

=  -  tt/5 

p  u  X. 

(III-38) 

n3 

=  -  sign(Vr)  p  Cd  £d(0) 

(III-39) 

55 


L 


■fa 


C1  ^  2  S1  3,2 


C2  4  S2  *2 


Xw  4  -  1/2  p  C£  u2  1 


(III-40) 


(III-41) 


(III-42) 


The  net  yaw  partial  derivatives  are  defined  by 


nij  variables. 


The  yaw  partials  with  respect  to  $  are 


3N1  3 

d<p  3<j>1 pw 


X  )  b 
sw  o 


(III-43) 


where 


JT  -  4  n2  V (0)  <p  (0) 


3T<Xpw  -  Xsw>  =  “  4  Xv 


The  yaw  partials  with  respect  to  p  are 


j—  s  JL(jj  4.  u  ) 
3p  3p 'Nsw  +  Nr 


(III-44) 


The  N  .  and  N  partial  derivatives  can  be  calculated  from 
sw  r 

equations  (11-44 (a))  and  (Il-ll(d)). 


The  yaw  moments  partial  derivatives  with  respect 
to  V  and  r  are  derived  similarly,  the  n^  coefficients 
related  to  these  modes  are 


9N 

119  9V 


(III-45) 


=  1M 

1111  “  9r 


(III-46) 


c.  Roll  Sensitivity  Coefficients 

The  roll  sensitivity  coefficients  corresponding 
to  state  matrix  elements  a^'s  are  determined  from  equation 
(III-7) .  They  are 

a71  =  k71/;Ck 

a73  *  k737lk 

a75  =  k76/Xk 

a76  =  k76/:Ek  +  n116/J 

a77  =  k79/Jk  +  n119/Z 

a711  =  k711/^Ik  +  nllll/:[ 


where  inertial  terms  I  and  1^  are  found  from  equations  (III- 
8(a))  and  (III-8(c)). 
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d.  Yaw  Sensitivity  Coefficients 

For  yaw  motion,  the  sensitivity  coefficients-  — 
aij's  have  been  determined  from  equation  (III-6) . 

These  coefficients  are  listed  below; 


alll  “  k71/2: 


a113  =  k73/;E 


a115  “  k75'^1 


a116  "  k76/3:  +  n116/^Ir 


a117  k79/;E  +  n119/^r 


allll  k711/^’  +  nllll/:Er 


where  the  inertial  term  Ir  is  calculated  from  equation 
(III-8 (b) ) . 


D.  STATE  SPACE  REPRESENTATION  AND  OUTPUT  EQUATION 

The  state  vector  of  the  linear  5  DOF  model  is  shown  as 


X  =  [dZ  dW  d0  dq  dM^  d<{>  dp  dy  dV  d^»  dr] 


Then,  the  system  state  equation  must  be 


X  =  Ax  +  Bu 


CIII-47) 


(III-48) 
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And  the  output  equation  is 


=  C  X 


(III-49) 


"A"  matrix  elements  have  already  been  determined.  However 
for  "B"  matrix  elements,  the  step  input  to  the  system  can 
be  considered  as  a  weight  removal  disturbance  and  for  the 
roll  and  yaw  modes,  a  rudder  order  disturbance. 

For  weight  removal  forcing  function,  the  input  is  to  be 


m 


(III-50) 


where 


Wf  *  final  craft  weight 
VT  =  initial  craft  weight 

m  =  craft  mass  corresponding  to  the  final 
weight  . 


When  the  rudder  order  is  used,  sway,  roll  and  yaw  forcing 
functions  are  generated  as  follows; 


sway 


roll 


b g i  =  i<2  xt  sin  6r  +  Ct  5r>  (III-51) 

b?1  =  (2  Xt  sin  5r  +  Cfc  6r>  (-  zr/Ik  +  xr/I) 


yaw 


bin  =  (2  Xt  sin  6r+Cfc  6r>  (xr/Ir  +  zr/I) 
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Finally,  defining  state  variables  x^,  etc., 
the  state  equations  of  the  5  DOF  linear  model  are  given 
as  follows;  ** 


X1  "  x2 


x2  ”  a21xl  +  a23x3  +  a24x4  +  a25x5 


x3  '  x4 


x4  =  a41xl  +  a43X3  *  a44x4  +  a45x5 


x5  *  a51xl  +  a55x5 


x6  "  x7 


x7  ’  a71xl  +  a73X3  +  a75XS  +  a76x6  +  a77x7 


+  a79x9  +  a7Uxll 


X8  =  x9 


x9  *  a96x6  +  a97x7  *  a99x9  *  a9UxU 


x10  -  XU 
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‘11 


alllXl  +  a113X3  +  a115x5  +  a116X6 


+  a117x?  +  aligx9  +  a1111x11 


All  coefficients  of  the  state  equations  are  defined 
in  the  preceding  sections. 

Considering  outputs  of  the  linear  model,  only  incremental 
changes  of  the  variables  can  be  observed. 
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IV.  LINEAR  FIVE  DEGREE-OF-FREEDOM 
^  -COMPUTER  PROGRAM  DESCRIPTIONS—  *  ' 

The  Digital  Simulation  Language  (DSL/360)  is  used  for 
the  linear  5  DOF  XR-3  computer  program.  Eleven  first  order 
integrator  blocks  are  needed  for  five  degree  of  freedom  simu¬ 
lations  and  for  all  linear  runs,  fourth  order  Runge-Kutta 
with  fixed  interval  integration  method  is  selected. 

In  addition  to  the  craft  dimensional  values,  initial  and 
final  steady  state  operating  points  are  entered  into  the 
program  for  linearization  purposes  and  plot  outputs. 

The  input  parameters  of  this  program  are 

V  =  Craft  speed  (ft/sec) 

WF  =  Final  craft  weight  (lb) 

W  ■  Initial  craft  weight  (lb) 

RUDANG  =  Rudder  input  (positive  for  port  turn) 

DR  =  Average  deadrise  angle 

PBBAR  =  Final  steady  gtate  value  of  plenum 
pressure  (lb/ft  ) 

ALD  =  Final  steady  state  craft  draft  (inches) 

PHI  »  Final  steady  state  roll  angle 

FYP  *  The  propulsion  system  thrust  corresponding 
to  the  craft  speed  (for  each  of  propellers) 

VL  »  Steady  state  lateral  speed  (ft/sec) 

R  *  Steady  state  yaw  rate  (rad/sec) . 
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The  steady  state  values  of  the  variables  are  determined 


from  6  DOF  nonlinear  model  simulations . 
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V.  SIX  DEGREE -OF -FREEDOM  NONLINEAR  MODEL ^SIMULATIONS 

In  this  chapter  6  DOF  nonlinear  model  simulations  are 
discussed. 

At  first, some  of  the  runs  are  made  in  order  to  determine 
craft  initial  steady  state  operating  point  corresponding  to 
craft  speed  and  thrust  inputs.  Since  6  DOF  model  weight 
transient  response  can  be  analyzed  by  means  of  entering  proper 
craft  initial  values  into  the  program,  this  is  simply  an 
initial  condition  response  of  the  nonlinear  dynamic  system. 

A.  WEIGHT  REMOVAL  TRANSIENT 

For  30  and  20  knot  craft  speeds,  the  steady  state  values 
of  some  variables  are  listed  in  Table  I. 


Table  I 

Steady  State  Values  of  the  Weight  Removal  Response 


Speed 

(knt) 

Pitch 

(degree) 

Pb 

(lb/ft^) 

Draft 

(inch) 

Thrust 

(lb) 

Weight 

(lb) 

30 

0.95 

24.8 

8.6 

355.6 

6722 

30 

0.26 

24.8 

5.6 

287.8 

6050 

20 

0.45 

24.8 

6.06 

186.2 

6050 

For  30  knot  craft  speed,  two  different  steady  state  pitch 
angles  are  observed  corresponding  to  the  different  weights. 

In  [2]  for  pitch  transient  response,  0.65°  steady  state  pitch 
angle  value  has  been  used  for  both  linear  and  nonlinear 
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simulations.  On  the  other  hand  as  seen  in  Table  I,  0.26° 
and  0.95°  steady  s£ate  pitcLfe,  vaiues  are  obtained.  This 'dis-  "* 
agreement  may  be  explained  by  noting  that  different  dimensional 
values  were  used  or  more  probably  a  different  modified  6 
DOF  computer  program  deck  was  used. 

For  6050  lb  craft  weight,  the  steady  state  pitch  angle 
0.26°  is  in  agreement  with  the  value  determined  by  Menzel  in 
figure  6  in  [6]. 

During  the  weight  removal  response,  simulations  of  6  DOF, 
two  different  stern  seal  base  leakage  areas  are  used.  In  [1] 
a  smaller  leakage  area  is  assumed  for  the  linear  heave  only 
model  as  A^  =  0.384.  However  A^  =  3.84  also  has  been  simulated 
and  it  is  understood  that  the  smaller  leakage  area  is  more 
reasonable  considering  other  steady  state  variables  in  the 
air  flow  dynamic  equations  (11-23)  and  (11-24) .  Since  at 
the  steady  state  condition  Qin  =  Qout  must  be  satisfied,  that 
permits  only  smaller  A^  for  a  given  Qq  and  steady  state 
values.  Therefore  A^  =  0.384  is  preferred  and  used  for 
nonlinear  simulations. 

B .  TURNING  MANEUVER 

When  the  3  DOF  pitch,  heave  and  roll  linear  model  was 
developed,  the  6  DOF  nonlinear  model  simulations  were  used 
in  order  to  determine  all  steady  state  roll  moments  during 
the  turn  maneuver.  In  addition  to  the  roll  moments,  yaw 
and  sway  outputs  were  also  investigated  and  the  results  are 
tabulated  in  Tables  II,  III  and  IV. 


Table  II 

Steady  State  Roll  Moments 
(rudder  =  5°,  30  knot  speed) 


Roll  Moments 
Stern  seal  FKSS 
Bow  seal  FKBS 
Sidewall  FKSW 
Rudder  FKRUD 
Propulsion  FKP 
Aerodynamic  FKAED 
Plenum  ABPB*PHI* (-Z) 


Positive 


131. 

0.057 

215.3 


Negative 

-2.2 

-280. 

-41.5 


-22.7 


Yaw  Moments 


Table  III 

Steady  State  Yaw  Moments 
(rudder  =  5°,  30  knot  speed) 

Positive 


Negative 


Stern  seal  FNBS 

3.22 

Bow  seal  FNSS 

0. 

0. 

Sidewall  FNSW 

- 

-482.2 

Rudder  FNRUD 

532.65 

- 

Propulsion  FNP 

2.78 

— 

Aerodynamic  FNAED 

- 

-  56.41 

Table  IV 


Steady  State  Sway  Forces  „ 
(rudder  =5°,  30  knot  speed) 

Sway  Forces  Positive 

Sidewall  FYSW 
Rudder  FYRUD 
Aerodynamic  FYAED 

Centrifugal  -R*U*AM  419.24 


Negative 

-351.8 

-  47.6 

-  19.07 


From  Table  II,  it  is  understood  that  for  a  given  (+) 
rudder  input,  the  sign  convention  of  the  craft  dynamics  was 
in  agreement  with  turning  to  the  left.  However,  a  negative 
sway  displacement  plot  is  observed  for  this  maneuver.  As  a 
result,  in  the  R.H.S.  subroutine  of  the  6  DOF  program,  the 
sign  correction  should  be  done  for  y-displacement . 

The  steady  state  sway  forces  acting  on  the  XR-3  test 
craft  in  a  turn  to  the  left  are  shown  in  figure  6.  As  seen 
in  this  figure  and  from  Table  IV,  the  rudder  sway  force  is 
in  the  steady  state  condition.  For  a  left  turn,  initially 
the  rudder  sway  force  is  generated  in  the  reverse  direction 
and  finally  a  negative  steady  state  rudder  sway  force  is 
produced.  Further  investigation  is  done  in  the  RUDDER  sub¬ 
routine  dynamic  equations  of  the  6  DOF  computer  model,  but 
no  error  has  been  determined  in  these  equations.  Also, 
different  craft  speeds  and  different  effective  rudder  angle 
coefficients  (ENDFAC)  are  used  but  everytime  the  same 
configuration  is  observed. 


In  addition  to  rudder  dynamics,  the  propulsion  system 
of  the  XR-3  contributes  to  the  roll,  yaw  and  sway  motions. 
However  this  contribution  is  not  a  dominant  factor  compared 
with  the  rudder  effect. 

For  the  XR-3  test  craft  propulsion  system  configuration, 
the  direction  of  the  thrust  vector  is  also  changed  as  well 
as  the  rudder  deflection.  As  seen  from  Tables  II,  III  and 
the  simulations  with  different  rudder  order,  it  is  understood 
that  the  propulsion  thrust  component  due  to  the  rudder  order 
is  not  calculated  properly,  therefore  the  following  modifica¬ 
tion  is  made  in  the  PROPULSION  subroutine  of  the  6  DOF  program, 

FYS  =  -STHSTS*CD+THSS*CD 

where  FYS  represents  the  sway  force  due  to  the  starboard 
propeller,  and  the  others  are  defined  in  the  propulsion 
subroutine  program. 

C.  DEADRISE  ANGLE  EFFECT 

As  reported  in  [4] ,  the  deadrise  angle  has  a  major  effect 
on  the  surface  effect  ship's  roll  behavior  and  therefore  is 
an  important  design  consideration.  The  deadrise  angle  effect 
on  the  roll  motion  has  been  investigated  by  Reidel  in  [3]. 

In  this  investigation,  some  arbitrary  coefficients  (PROMO  1, 

2,  etc.)  are  used  for  the  vertical  force  component  of  the 
sidewall  hydrodynamics.  From  this  investigation,  the  conclusion 
reported  by  Reidel  was  that  more  damping  could  be  obtained  for 
the  pitch  and  roll  motions  by  changing  the  deadrise  angle. 
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The  deadrise  effect  is  reinvestigated  by  simulating  the 
system  with  different  rudder  orders  including  zero  disturbance, 
then  the  heave  forces  due  to  sidewall  hydrodynamics  were 
compared  to  each  otheri.  As  a  result,  no  significant  sidewall 
heave  force  change  was  observed.  For  instance,  for  0°  and 
5°  rudder  inputs,  the  sidewall  heave  forces  (FZSW)  are  found 
to  be  -568.8  and  -568.05  respectively.  This  result  is  not 
acceptable,  since  in  a  turn  maneuver,  the  cross  flow  drag 
is  generated  at  the  outboard  sidewall,  so  the  deadrise  projec¬ 
tion  of  cross  flow  drag  must  be  contributed  to  the  sidewall 
heave  force.  More  investigation  has  been  made  in  the  SIDEWALL 
subroutine  of  the  6  DOF  program  and  it  has  been  determined 
that  only  the  force  component  due  to  added  mass  given  by 
equation  (11-21)  is  projected  as  a  lift  force.  Actually  the 
cross  flow  drag  component  must  be  included  to  this  projection, 
since  the  frontal  area  for  cross  flow  drag  component  is  out¬ 
board  of  the,  submerge.d  section  o.f  the  sidewall.  Therefore 
it  is  concluded  that  the  XR-3  test  craft  6  DOF  program  SIDEWALL 
subroutine  does  not  include  effectively  the  lift  force  due 
to  deadrise.  On  the  other  hand,  the  lift  force  and  horizontal 
drag  component  can  have  different  drag  coefficients,  therefore 
the  modification  of  the  sidewall  subroutine  is  left  for  future 


studies. 


VI.  COMPARISON  OF  NONLINEAR  AND  LINEAR  MODEL  RESULTS 

A.  C.G.  ACCELERATION  (see  figs.  9  and  23) 

A  10%  weight  removal  disturbance  is  used  for  both  models. 

The  linear  and  nonlinear  model  showed  very  rapid  response 
with  small  rise  time.  As  seen  in  figure  9,  the  6  DOF  model 
transient  is  relatively  more  damped  compared  with  the  linear 
model  response.  The  linear  model  showed  a  small  overshoot 
but  it  is  also  sufficiently  damped.  Small  oscillatory  behavior 
is  observed  in  both  figures  which  is  due  to  pitch-heave 
coupling.  In  general  the  two  models  agree  in  transient  response 
and  steady  state  behavior. 

Both  of  the  C.G.  acceleration  curves  were  normalized  by  a 
factor  “g'^the  gravitational  constant. 

B.  PITCH  TRANSIENT  (see  figs.  10  and  24) 

As  seen  in  these  figures,  both  models  showed  an  oscillatory 
transient  response  with  the  same  frequency,  w  =  3.12  rad/sec. 
However  their  steady  state  values  are  different.  In  fact, 
the  pitch  moment  due  to  plenum  pressure  and  the  effective 
plenum  pressure  center  play  a  dominant  role  for  pitch  motion. 

In  [2]  Barnes  developed  plenum  modeling  including  effective 
plenum  roof  area  calculations  and  plenum  pressure  centroid 
shift  due  to  draft  variation.  Therefore  the  linear  model 
plenum  modeling  is  different  from  the  nonlinear  modeling, 
so  that  a  greater  plenum  moment  arm  is  calculated  in  the  linear 
model. 
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The  steady  state  pitch  angle  values  are  0.5°  and  0.75° 
for  nonlinear  and  linear  models  correspondingly. 

C.  PLENUM  PRESSURE  TRANSIENT  (see  figs.  7  and  21) 

As  seen  in  these  figures,  both  models  are  exactly  in 

agreement  with  each  other.  In  each  curve,  the  plenum  pressure 

2 

suddenly  drops  to  21.6  (lb/ft  )  then  increases  up  to  its 
steady  state  value  with  sufficient  damping.  This  is  expected 
behavior  because  the  plenum  volume  increases  when  10%  weight 
is  removed,  obviously  this  disturbance  causes  a  sudden  drop  in 
plenum  pressure.  As  seen  before,  pitch-plenum  pressure  coupling 
had  been  observed  during  the  exponential  increasing  phase  of 
the  plenum  pressure.  At  the  fifth  second  of  simulation  time, 
nonlinear  model  plenum  pressure  is  about  24.29  while  the 
linear  model  pressure  is  24.25. 

D.  DRAFT  TRANSIENT  (see  figs.  8  and  22) 

Excellent  “agreement  was  achieved  between  linear  and  non¬ 
linear  model  draft  responses.  As  expected,  the  draft  decreases 
exponentially  when  the  weight  is  removed  from  the  craft. 

The  draft  is  found  to  be  5.7  inches  for  the  linear  model, 

5.85  inches  for  the  nonlinear  model  at  the  fifth  second  of 
simulation  time. 

As  a  result,  both  models  exhibited  the  same  draft 
responses. 

E.  ROLL  RESPONSE 

As  mentioned  before,  at  first  a  3  DOF  linear  model  (pitch- 
heave-roll)  had  been  developed.  3  DOF  linear  and  6  DOF 
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nonlinear  models  roll  behavior  comparison  is  done  by  means 
of  a  step  input  roll  moment  disturbance  to  the  linear  model 
and  a  turning  command  to  the  nonlinear  model.  In  figure  13 
6  DOF  nonlinear  model  roll  angle  response  is  shown  for  a  5° 
rudder  order  (left  turn )/  in  this  simulation  1.18°  outboard 
steady  state  roll  angle  is  generated.  From  the  R.H.S.  sub¬ 
routine  of  the  6  DOF  program,  all  the  steady  state  roll  moment 
components  are  determined  as  seen  in  Table  II  and  the  stun  of 
these  disturbance  moments  applied  to  the  linear  model  roll 
mode  as  a  step  input.  Roll  behavior  is  analyzed  from  figure 
27.  As  expected,  the  linear  model  showed  a  greater  overshoot 
and  more  oscillatory  response  compared  with  the  nonlinear 
roll  transient.  Since  the  linear  model  has  been  disturbed 
by  a  step  roll  moment  corresponding  to  the  nonlinear  model 
steady  state  roll  moments.  On  the  other  hand,  its  steady 
state  roll  angle  is  found  to  be  1.28°.  The  roll  steady  state 
difference  of  two  models  is  about  8.4%.  For  a  3°  rudder  input, 
nonlinear  and  linear  simulations  are  obtained  in  figures  11 
and  25  with  the  steady  state  values  0.62°  and  0.65°  respec¬ 
tively.  Hence  the  difference  is  4.8%  in  the  steady  state 
values  for  these  simulations. 

Figures  12  and  14  show  XR-3  yaw  rate  plots  and  corres¬ 
ponding  linear  simulations  are  shown  in  figures  26  and  28. 

As  a  result,  the  3  DOF  linear  model  and  6  DOF  program 
roll  simulations  are  in  agreement  with  each  other  within 
10%  difference. 
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As  discussed  in  the  previous  chapter,  some  problems 
are  observed  in  the  6  DOF  nonlinear  program.  Therefore 
developing  two  more  modes  for  the  3  DOF  model,  the  5  DOF 
linear  model  has  been  generated.  By  means  of  this  model, 
roll  behavior  of  the  XR-3  in  a  turn,  deadrise  ingle  and 
cross  product  inertial  term  effects  have  been  investigated 
for  20  knot  craft  speed  with  different  rudder  deflections. 

The  roll  response  can  be  seen  in  figure  29  omitting  the 
cross  product  inertial  term  (I  =  -28Q0.)  in  the  equations 
of  motion,  applying  a  10°  rudder  order  and  assuming  70  degree 
average  deadrise  angle  of  the  craft.  For  this  particular 
test,  the  craft  starts  inboard  rolling  with  the  oscillatory 
transient  then  goes  to  0.47°  steady  state  roll.  Using  the 
same  input  parameters  and  also  including  the  cross  product 
term  I  in  the  equations  of  motion,  smaller  overshoot  (18%) 
and  a  more  damped  roll  response  has  been  observed  in  figure 
31..  £s  seen  in  this  figure,  at  the  initial  phase,  secondary 
overshoot  has  been  generated.  The  steady  state  roll  angle  is 
0.40°  and  the  oscillation  frequency  w  =  3.5  rad/sec.,  for 
this  simulation.  On  the  other  hand,  the  corresponding  6  DOF 
nonlinear  simulation  is  shown  in  figure  17  and  also  a  roll 
rate  plot  is  shown  in  figure  18.  The  nonlinear  simulation 
in  figure  17  shows  a  more  damped  transient  but  greater  steady 
state  roll  behavior  (1.32°)  compared  with  the  corresponding 
linear  simulation  shown  in  figure  31. 
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Figures  30  and  32  present  linear  simulation  curves  for 
5  and  15  degree  rudder  inputs  respectively.  Both  curves 
exhibit  the  same  linear  model  response  discussed  before. 
Finally  the  linear  and  nonlinear  model  steady  state  responses 
are  summarized  as  follows. 


Rudder 

6  DOF 
Nonlinear 

5  DOF 
Linear 

Reference 

5 

0.57 

0.24 

figures  15  and  30 

10 

1.32 

0.40 

figures  17  and  31 

15 

2.2 

0.45 

figures  19  and  32 

The  deadrise  angle  effect  is  investigated  by  using  differ¬ 
ent  deadrise  angle  values  for  10  degree  rudder  input  and  20 
knot  craft  speed  simulations.  The  results  are  tabulated  in 
Table  V. 
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For  the  linear  model,  the  lateral  velocity  and  yaw  rate 
plots  are  shown  in  figures  33  and  34.  Both  of  them  showed 
sufficiently  damped  transient  response.  From  figure  33, 
the  lateral  speed  steady  state  value  is  found  3.6  (ft/sec) 
and  as  expected  for  turning  to  the  left,  the  negative  yaw 
rate  is  generated  as  -0.1  (rad/sec). 
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Figure  7.  6  DOF  Plenum  pressure  transient 
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Figure  8.  6  DOF  draft  transient 
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Figure  9.  6  DOF  C.G.  Acceleration 
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Figure  10.  6  DOF  pitch  angle  transient 
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Figure  11.  6  DOF  Roll  angle  plot,  5  =  30° 
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Figure  12.  6  DOF  roll  rate  plot,  6  =  3< 
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Figure  13.  6  DOF  roll  angle  plot,  6  =  5° 
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Figure  14.  6  DOF  roll  rate  plot,  &  =  5° 
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Figure  15.  6  DOF  roll  angle  plot,  6  =5° 
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Figure  16.  6  DOF  roll  rate  plot,  6 
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Figure  17.  6  DOF  roll  angle  plot,  6r 
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Figure  18.  60  DOF  roll  rate  plot,  5  =  10° 
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Figure  19.  6  DOF  roll  angle  plot,  6 
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Figure  20.  6  DOF  roll  rate  plot,  5  =  15° 
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Figure  21.  Linear  5  DOF  plenum  pressure  transient 
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Figure  22.  Linear  5  DOF  draft  transient 


Figure  23.  Linear  5  DOF  C.G.  acceleration 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  the  preceding  sections/  the  5  DOF  linear  XR-3  model 
was  developed/  then  weight  removal  response  and  roll  behavior 
are  compared  with  6  DOF  nonlinear  simulations.  Also  some 
of  the  craft  sensitivity  parameters  are  investigated  for  the 
XR-3  test  craft. 

The  results  of  this  study  are  summarized  and  recommenda¬ 
tions  are  made  as  follows. 

1.  The  system  response  is  investigated  by  both  a  linear 
model  and  a  6  DOF  nonlinear  model  by  means  of  sudden  weight 
removal  from  the  C.G.  and  by  sway  force  disturbance  in  a 
turn. 

Comparison  of  results  shows  that  the  difference  is 
within  10%  for  both  models.  In  fact,  assuming  accuracy  in 
measurements  of  about  10%,  the  linearized  model  gives 
acceptable  results  as  compared  with  the  nonlinear  model. 
Obviously  the  linear  model  simulations  take  less  computer 
CPU  time. 

2.  Concerning  roll  response  during  the  turn  maneuver, 

the  product  inertia  term  (I  )  affects  transient  response 

xz 

such  as  overshoot,  damping,  etc.  More  overshoot  and 
oscillatory  response  have  been  observed  without  Ixz* 

However,  this  term  is  found  to  be  -2800.  from  the  6  DOF 
program.  It  is  suggested  that  this  number  be  verified. 
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3.  The  importance  of  deadrise  is  observed  for  roll 
steady  state  behavior.  The  XR-3  test  craft  steady  state 
roll  angle  depends  upon  the  average  deadrise  angle  of  the 
sidewall.  Different  average  deadrise  angles  have  been 
simulated  by  means  of  the  linear  model  and  it  is  observed 
that  the  craft  rolls  inboard  or  outboard  according  to  the 
sidewall  deadrise  configuration. 

The  6  DOF  nonlinear  model  SIDEWALL  SUBROUTINE 
must  be  modified  by  including  a  cross  flow  drag  component 
deadrise  projection  (lift  force) .  However  a  new  drag 
coefficient  is  needed  for  this  lift  component. 

4.  The  aerodynamic  effect  of  the  craft  contributes 
about  5%  of  the  total  yaw  moment. 

5.  The  XR-3  failure  mode  operation  (such  as  fan  failure) 
can  be  simulated  by  using  either  the  linear  or  nonlinear 
model.  But,  the  nonlinear  model  is  preferable  for  failure 
mode  simulations,  since  it  includes  sidewall  and  stern 
seal  leakage  configurations. 
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